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^ . ABSTRACT 



We compare the semi-analytic models of galaxy formation of Fu et al. (2010), which 
track the evolution of the radial profiles of atomic and molecular gas in galaxies, with 
gas fraction scaling relations derived from a stellar mass-limited sample of 299 galaxies 
O " from the COLD GASS survey. These galaxies have measurements of the CO(l-O) line 

from the IRAM 30-m telescope and the HI line from Arecibo, as well as measurements 
^ ■ of star formation rates, stellar masses, galaxy sizes and concentration parameters from 

^ " GALEX+SDSS photometry. The models provides a good description of how condensed 

baryons in star-forming galaxies are partitioned into atomic and molecular gas and 
stars as a function of galaxy stellar mass and stellar surface density. The models do not 
J> I reproduce the observed tight relation between stellar surface mass density and bulge- 

. to-disk ratio in these galaxies. The current implementation of "radio-mode feedback" 

in the models produces trends that disagree strongly with the data. In the models, 
gas cooling shuts down in nearly all galaxies in dark matter halos above a mass of 
~ IO-'^^Mq. As a result, stellar mass is the observable that best predicts whether a 
" galaxy has little or no neutral gas, i.e. whether a galaxy has been quenched. In contrast, 

, our data show that quenching is largely independent of stellar mass. Instead, there are 

clear thresholds in bulge-to-disk ratio and in stellar surface density that demarcate 
the location of quenched galaxies in our chosen parameter space. We speculate that 
the disagreement between the models and the observations may be resolved if radial 
transport of gas from the outer disk is included as an additional bulge-formation 
KnJ I mechanism in the models. In addition, we propose that processes associated with 

^ ■ bulge formation play a key role in depleting the neutral gas in galaxies and that gas 

accretion is suppressed in a significant fraction of galaxies following the formation of 
the bulge, even in dark matter halos of low mass. 
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- radio lines: galaxies - surveys 



1 INTRODUCTION 

An important goal in modern galaxy formation is an im- 
proved understanding of the physical processes that regulate 
the rate at which stars form in galaxies. These processes in- 
clude the cooling and accretion of gas within dark matter 
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halos, the transformation of the accreted gas into molecular 
clouds and stars, and the effects of "feedback" from massive 
stars and accreting black holes on the gas in and around 
galaxies. 

The Lambda-Cold Dark Matter model provides de- 
tailed predictions for how dark matter halos assemble over 
time. Within these halos, gas will cool and settle into a 
rotationally-supported disk. One hotly debated issue is the 
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degree to which gas loses angular momentum to the sur- 
rounding dark matter during this process. In many "semi- 
analytic" models, the angular momentum of the gas is as- 
sumed to bo conserved. This produces a size-circular velocity 
relation for galactic disks that is in good agreement with ob- 
servations (Mo ,Mao & White 1998). More detailed numeri- 
cal hydrodynamical simulations of disk galaxy formation in 
a ACDM universe demonstrate, however, that the amount 
of angular momentum that is lost during the collapse of 
the dark matter halo depends sensitively on how feedback 
processes are included in the models (Navarro & Steinmetz 
1997; Governato et al 2004,2007). In many simulations, the 
disks form early and are too compact. The resulting star for- 
mation timescales are then too short to be consistent with 
observations (e.g. Oser et al 2010). 

Another important unsolved problem is to understand 
why galaxies divide into two distinct "families" - those with 
ongoing active star formation, and those where star forma- 
tion has been quenched (Strateva et al 2001; Kauffmann ct 
al 2003b; Baldry et al 2004). It has now become clear that 
these distinctions existed already at redshifts as high as 2.5 
(Williams ct al 2010; Wuyts ct al 2011), when cold gas ac- 
cretion rates in almost all dark matter halos are predicted 
to be very high. 

Many recent theoretical gala:xy formation models that 
aim to reproduce the statistical properties of the massive 
galaxy population (e.g. Croton ct al 2006; Bower ct al 2006; 
Cattaneo et al 2006 De Lucia & Blaizot 2007; Somerville et 
a 2008; Guo et al 2011; Lu et al 2011) invoke star-formation 
quenching mechanisms that set in at a characteristic dark 
matter fialo mass. The characteristic mass is associated with 
the transition between the regime where gas cooling times 
are short compared to the dynamical time of the dark matter 
halo and gas accretes in the form cold, condensed clouds, and 
the regime where cooling times are long and gas accretes 
from a corona of gas that is in virial equilibrium with the 
surrounding halo. 

One might ask why lowered rates of star formation in 
galaxies should be linked with this transition. One argument 
that is often used is that rclativistic jets, generated when gas 
accretes onto black holes, heat the surrounding hot gas and 
prevent it from forming stars (see, for example, Croton et 
al. 2006). The conditions under which black holes produce 
jets are not well understood. There is clear observational ev- 
idence that radio-emitting jets play a role in regulating the 
cooling of gas in nearby gala^xy clusters (see McNamara & 
Nulsen 2007 for a recent review). Theoreticians have thus 
made an "ansatz" that jets from radio-loud AGN will pre- 
vent gas from cooling and forming stars in all dark matter 
halos that are predicted to contain a hot gas atmosphere. 
The characteristic halo mass threshold that separates star- 
forming and quiescent galaxies is the main factor that deter- 
mines how massive galaxies evolve in current models, so it 
is clearly important to test the existence of such a threshold 
using real observations. 

We note that accretion and quenching processes affect 
the gas components of galaxies. The most direct empirical 
constraints on how they operate thus come from observa- 
tions of the gas in galaxies. This has been the primary mo- 
tivation for the Galex Arecibo SDSS survey (GASS), as well 
as the CO Legacy Database for GASS (COLD GASS), which 
are measuring the atomic and molecular gas contents of an 



unbiased sample of several hundred galaxies with redshifts 
between 0.025 and 0.05, and with stellar masses in the range 
10^° < M» < 1O"-^M0. In both surveys, the strategy is to 
observe each galaxy until the HI/CO lines are detected, or 
until upper limits in the ratio of atomic and molecular gas 
mass to stellar mass of ~ 0.015 reached. The aim of the 
program is to carry out a census of the condensed baryons 
in galaxies in the local Universe, to study scaling relations 
between the gas and stellar properties of these galaxies, and 
to understand gas accretion and quenching in galajcies. 

Details of the GASS and COLD GASS survey designs, 
as well as target selection and observing procedures are 
given in Catinella ct al (2010) and Saintonge ct al (2011a). 
These two papers also presented relations between the HI 
and H2 mass fractions of galaxies (defined as Mhi/M, and 
Mhj/M*) and global galaxy parameters such as stellar mass 
M*, stellar mass surface density /i,, galaxy bulge-to-disc ra- 
tio (as parameterized by the concentration index C of the 
r-band light), and specific star formation rate SFR/M« (see 
also Schiminovich et al 2010). The two surveys uncovered 
sharp thresholds in /u* and C below which all galaxies have a 
measurable cold gas component, but above which the detec- 
tion rate of the CO and HI lines drops suddenly, suggesting 
that "quenching processes" have occurred in these systems. 

In this paper, we compare the observed gas fraction re- 
lations from the COLD GASS survey with predictions from 
the semi-analytic models of Fu ct al (2010; hereafter FIO), 
which track the formation of molecular gas in disk galajcies. 
The FIO models are based on the L-galaxies semi-analytic 
code, described in detail in Croton et al (2006) and updated 
in Dc Lucia & Blaizot (2007). 

The models are currently being configured to operate 
on the latest update of the L-galaxies code by Guo et al. 
(2011); the comparison in this paper will be restricted to 
the version of the models in the published FIO paper. The 
main new aspect of these models is that galactic discs are 
represented by a series of concentric rings in order to track 
the evolution in the gas and stellar surface density profiles 
of galaxies over cosmic time. Two simple prescriptions for 
molecular gas formation processes are included: one is based 
on the analytic calculations by Krumholz, McKee & Tum- 
linson (2009; hereafter KMT), and the other is a prescrip- 
tion where the H2 fraction is determined by the pressure of 
the interstellar medium (Blitz & Rosolowsky 2006, hereafter 
BR), with an implementation similar to that in Obrcschkow 
et al. (2009). The free parameters of the models that regulate 
the rate at which gas is turned into stars and the efficiency 
with which supernovae reheat g function of halo mass 

have been tuned to reproduce a number of key observables, 
including the observed present-day galaxy luminosity func- 
tion, the gas-phase mass-metallicity relation, and the mean 
HI gas fraction as a function of B-band luminosity (see De 
Lucia, Kauffmann & White (2004) for a more detailed dis- 
cussion) . 

There have been a number of other cosiiiological mod- 
els that follow the formation of molecular gas in disks. La- 
gos et al (2011) have developed semi-analytic models that 
track molecular gas formation in cosmological simulations 
of galaxy formation, but these models do not predict the 
detailed radial profiles of galaxies as in FIO. Robertson & 
Kravtsov (2008) and Dutton et al (2010) model the radial 
profiles of the gas and stars in galaxies, but their models 
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are not embedded within full cosmological N-body simula- 
tions. Gnedin et al (2009) have carried out very high reso- 
lution cosmological simulations of molecular gas formation 
that include radiative transfer and detailed treatment of the 
chemistry, but their box sizes are too small to make statis- 
tical predictions. 

The FIG paper demonstrated that their models could 
fit the radial HI, H2, stellar mass and SFR profiles of spi- 
ral galaxies from the THINGS/HERACLES surveys. In this 
paper, we investigate whether the models also yield correct 
gas fraction scaling relations and distribution functions. As 
will be seen, it is useful to break this analysis into two dis- 
tinct parts: l)comparison between models and observations 
of the distribution of gas fractions in the population of galax- 
ies with detectable gas, 2) comparison of how absence of gas, 
which we refer to as "quenching" , is manifested as a function 
of observable quantities such as stellar mass, stellar surface 
density, and bulge-to-disk ratio. As we will show, such a com- 
parison elucidates those aspects of the models which work 
well, and those that fail. Throughout this paper, we have as- 
sumed a cosmology with Q, — 0.3, A — 0.7 and Hq — 70 km 
s~^ Mpc~^ to derive our observational quantities. The cos- 
mology assumed for the FIO model is a ACDM models with 
n = 0.25, A = 0.75, Ho = 73 km s"^ Mpc"^ and as = 0.9. 



2 GENERATING COMPARISON SAMPLES 
FROM THE SIMULATIONS 

The semi-analytic models used in this paper are nearly iden- 
tical to those described in detail in FIO. We have made only 
one small change. The free parameters in the FIO models 
were chosen so that reasonable fits to a variety of observ- 
ables could be obtained. These observables included the H2 
mass function of Keres et al (2003) derived using galaxies 
drawn from the FCRAO Extragalactic CO Survey (Young et 
al 1995). These authors adopted a CO-to-H2 conversion fac- 
tor that was a factor of 1.5 times larger than the one adopted 
by Saintonge et al (2011a). To bring the simulations back 
into agreement with the H2 data using the conversion factor 
adopted by Saintonge et al, we simply increase the super- 
nova reheating rate (given by the parameter edisc) in Table 
1 of FIO by a similar factor (e^isc = 5 instead of 3.5). We 
note that changing edisc simply shifts the amplitude of the 
gas mass fraction scaling relations up or down. It does not 
change the predicted slope or scatter in these relations. We 
also note that the H2 mass function derived from the COLD 
GASS data for galaxies with logM* > 10^° Mq agrees well 
with that of Keres et al (2003) at the high mass end, once 
the same CO-to-H2 conversion factor is adopted (Fu et al 
2012, in preparation). 

The observational sample is an update to that described 
in Saintonge et al 2011(a), consisting of 299 galaxies with 
CO(l-O) line observations with the TRAM 30 m telescope. 
HI line measurements from the GASS survey (Catinella et al 
2010) are available for 270 out of the 299 galaxies. Because 
the observational sample is selected only by stellar mass, 
it is easy to compare simulation results directly with the 

^ c«CO = 3.2M0 (K km pc'^)"^, which does not include a 
correction for the presence of helium. 



data. There are only two observational selection issues that 
require careful treatment: 

(i) In both the GASS and COLD GASS surveys, targets 
are selected so that the resulting stellar mass distribution is 
approximately flat. This was done in order to be able to com- 
pute gas fractions in different stellar mass bins with roughly 
the same errors. The sampling must be taken into account 
when we look at gas fraction distributions as a function of 
parameters other than stellar mass. In this paper, we adopt 
two approaches: a) We either weight each galaxy in the sam- 
ple by the inverse of /(Af, ), where /(M«) is the "sampling 
rate function" required to transform the true stellar mass 
distribution to a flat one, b) We apply the sampling rate 
function /(M*) to the simulated galaxies to create "mock 
catalogue" with flat stellar mass distributions that can be 
compared directly with the data. 

(ii) It is important to take account of the detection lim- 
its of the survey when comparing model galaxies with the 
data. This is illustrated in the top panels of Figure 1, where 
we plot M(HI)/M, (left) and M(H2)/M, (right) as a func- 
tion of stellar mass for 299 galaxies that have been observed 
in CO as of July 2011 (270 of these have HI observations 
from the GASS survey). Galaxies where the HI or CO line 
was detected (see Catinella et al 2010 and Saintonge et al 
2011a for details about line detection procedures) are plot- 
ted in blue. Galaxies without detections are plotted in red 
at the positions of their 5a upper limits in M(HI)/M* or 
M(H2)/A'/.. 

To a very close approximation, the HI mass fraction limit 
log [M(HI)/M,]Hm is -1.82 (corresponding to a HI mass 
fraction limit 0.015 ) for galaxies with log A/* > 10.3. For 
galaxies with 10 < log A/* < 10.3, log [M(HI)/M*],i„ = 
-1.066 log Af. +9.16. 

The CO line detection limits are somewhat more com- 
plicated. As discussed in Saintonge et al (2011a), the inte- 
gration times are set so that log [M(H2)/M,]Km=-1.82 for 
galaxies with logAf* > 10.6. For lower mass galaxies, in- 
tegration times are nominally set to reach a fixed r.m.s. 
of around 1.1 mK per 20 km/s wide channel, but the ac- 
tual value fiuctuated according to weather conditions at 
the time of observation. Therefore, when creating our mock 
catalogues from the simulations, we simply impose a dis- 
tribution of upper limits similar to that seen in the top 
right panel of Figure 1. We take [M{R2) /Mt]iim.,max, the 
maximum possible value of the H2 mass fraction limit, as 
log [M(H2)/M,],,„,„,, = -1.72 if logAf. > 10.6, and log 
[MiUi) /M^]iim,ma^ = 8.78- log M, if log Af, < 10.6. The Ha 
mass fraction limit that we adopt is randomly distributed 
between the maximum value and a minimum value that is 
0.35 dex smaller. 

In this way, we are able to classify galaxies in the simu- 
lations as HI/CO "detections" or "non-detections" and com- 
pare properties such as stellar masses, stellar surface densi- 
ties, specific star formation rates and concentration indices 
with those of the real galaxies in our survey. We note that 
stellar masses and star formation rates are standard outputs 
of semi-analytic models. Because the FIO models are able to 
track radial profiles of the gas and the stars in galaxies, they 
are also able to predict stellar surface mass densities, which 
are defined as 0.5Af»/(7r7?§o)i where R50 is the radius con- 
taining half the stellar mass of the galaxy. 
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Figure 1. In the left panels, HI mass fraction (log M(HI)/M*) 
and specific star formation rate (log SFR/M,) are plotted as a 
function of the log M* . Blue points denote galaxies with HI line 
detections. Red points denote galaxies where HI was not detected. 
In the right panels, H2 mass fraction (log M(H2)/M*) and specific 
star formation rate (log SFR/M*) are plotted as a function of the 
log M« . Blue points denote galaxies with CO line detections. Red 
points denote galaxies where CO was not detected. 



In semi-analytic models, bulges form either as a result 
of a merger between two galaxies, or as a consequence of 
instabilities that set in when disks reach a certain critical 
threshold density. In the models, the total mass of the bulge 
is more reliably predicted than its size. As shown in Fig- 
ure 1 of Weinmann et al (2009), there is a tight correlation 
between the concentration index C (defined as the ratio of 
the radii containing 90% and 50% of the r-band light) and 
galaxy bulge-to-disk radio derived using 2-dimensional de- 
composition codes (Gadotti 2009). We have therefore elected 
to transform the bulge-to-disk ratios predicted by the mod- 
els to concentration indices using a fit to this relation: 
C=1.8+2.UB/T. 



3 RESULTS 

In the analysis presented in this paper, we will attempt to 
answer two questions: 

(i) Can the simple disk formation models described in the 
FIO paper explain observed gas fraction scaling relations for 
galaxies with gas? 

(ii) Does the transition between the population of galax- 
ies with gas and the population without gas occur in the 
same way in the models and in the data? 

We have chosen the gas fraction detection threshold of 
our two surveys as the nominal division between the popu- 
lation of galaxies we will henceforth refer to as "active" and 
the population that we will call "quenched" . Note that it 
is possible that some fraction of galaxies classified as "ac- 
tive" are actually losing their gas and evolving into quenched 



systems. Conversely, some fraction of galaxies classified as 
"quenched" may have simply run out of gas just prior to the 
next accretion episode. Such galaxies may be thought of as 
"transition" systems. 

In the bottom panels of Figure 1, we plot specific star 
formation rate versus stellar mass for the galaxies in our 
sample. In the left panel, galaxies with HI detections are 
colour-coded blue and those where the HI line was not de- 
tected are colour-coded red. In the right panel, we do the 
same thing according to whether the CO line was detected 
or not. As can be seen, galaxies with log SFR/M, < —11 
are not usually detected in HI and almost never in CO. 
This means that there is a rather clean distinction between 
"active" and "quenched" galaxies that applies to both gas 
fraction and specific star formation rate for the majority of 
galaxies in our sample. In section 3.2.2, we return to the is- 
sue of whether we can identify a sub-population of galaxies 
that is in the process of being quenched. 

3.1 Gas scaling relations for the active population 

In Figures 2 and 3 we compare the mean HI and H2 gas frac- 
tions as a function of M, , p*, C and SFR/M, for data and 
for models. The reader is referred to Catinella et al (2010) 
and Saintonge et al (2011a) for details on how HI and CO 
line fluxes are transformed into HI and H2 masses. We note 
that we have adopted a constant Galactic conversion factor 
ctco = 3.2Mq (K km s~^ pc^)~^, which does not include a 
correction for the presence of helium. A constant conversion 
factor should be a reasonable approximation for the galax- 
ies in our survey, which all have stellar masses greater than 
10^" M0 and gas- phase metallicities near solar (Saintonge et 
al 2011b). 

In Figures 2 and 3, the blue points show the mean values 
of M(HI)/M, and M(H2)/M, calculated using the galaxies 
with HI and CO line detections. Error bars are calculated 
using bootstrap resampling. The red and black curves in- 
dicate the mean values calculated from the model galaxies 
that are classified as detections (see section 2). The black 
curves show results from the FIG model that uses the KMT 
prescription for the conversion of atomic to molecular gas. 
The red curves show results from the FIO model that uses 
the BR pressure-based prescription. We note that there are 
many more simulated galaxies than observed galaxies. In 
generating the curves, we have ordered the galaxies by in- 
creasing M, C and SFR/M* and computed means for 
bins of 70 galaxies. 

Figures 2 and 3 show that the models reproduce the 
observed scalings between mean HI/H2 mass fraction and 
galaxy stellar mass, surface mass density and concentra- 
tion quite well. The differences between the HI mass frac- 
tions predicted by the KMT and BR prescriptions are small. 
Given the systematic uncertainties in the calibration of the 
BR prescription, the fact that the KMT prescription is 
based on simple analytic calculations, the fact that our semi- 
analytic disc formation models are extremely idealized, and 
the fact that the model free parameters were not tuned to 
fit the COLD GASS scaling relations, we find it remarkable 
how well these results agree with the data. 

The relation between H2 mass fraction and specific star 
formation is shallower than predicted by the models, which 
lie on a relation with slope unity. The discrepancy in the re- 
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Figure 2. The mean HI mass fraction (log M(HI)/M,) for COLD 
GASS galaxies with HI line detections (blue points) is plotted as 
a function of stellar mass, stellar mass surface density, concentra- 
tion and specific star formation rate. Error bars are calculated us- 
ing bootstrap resampling. Results for the FIO models are shown as 
black (KMT atomic-to-molecular gas prescription) and red (BR 
pressure prescription) curves. 
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Figure 3. The mean H2 mass fraction (log M(H2)/M») for 
COLD GASS galaxies with CO line detections (blue points) is 
plotted as a function of stellar mass, stellar mass surface density, 
concentration and specific star formation rate. Error bars are cal- 
culated using bootstrap resampling. Results for the FIO models 
are shown as black (KMT atomic-to-molecular gas prescription) 
and red (BR pressure prescription) curves. 



lation between M(H2)/M« and SFR/M» is easily understood 
in light of the results presented in Saintonge et al (2011b). 
The FIO models adopt the assumption that Espr — aSn, 
with a constant (corresponding to an effective molecular gas 
depletion timescale of 2 Gyr) in all galaxies. This assumption 
was motivated by the results in Leroy et al (2008). Saintonge 
et al (2011) showed, however, that there was more than a 
factor of 5 variation in the global molecular gas depletion 
timescale in different galaxies, and that the least actively 
star-forming galaxies (i.e lowest values of SFR/A/*) had the 
longest depletion times. This explains the shallow observed 
relation in the bottom right panel of Figure 3. 

In Figure 4, we show scatter plots of log SFR/M, ver- 
sus log M(H2)/M, and log M(HI)/M, for our survey galax- 
ies (blue points) and for galaxies from our mock catalogue 
(black points). Results are only shown for the KMT atomic- 
to-molecular gas prescription, because results for the pres- 
sure prescription are virtually the same. As can be seen the 
models reproduce the scatter in specific star formation rate 
versus atomic gas fraction quite well, but not in specific star 
formation rate versus molecular gas fraction. Clearly, there 
are additional processes at work that determine the rate at 
which molecular gas will form stars. This is the subject of a 
future paper (Saintonge et al, in preparation). 

As well as looking at mean HI and H2 mass fraction 
relations for galaxies with detections, it is instructive to an- 
alyze full distribution functions of gas mass fractions. These 
are illustrated in Figures 5 and 6. The data are shown as blue 
histograms, with error bars computed from bootstrap resam- 
pling. Black and red curves are for the models, as before. We 
plot distributions of log M(HI)/M, and log M(H2)/M. in 
three different intervals of stellar mass, stellar mass surface 
density and concentration. These plots show that the models 



provide a good representation not only of the mean HI/H2 
mass fractions, but also of the scatter around the mean. As 
discussed in FIO and will be illustrated below, there is a tight 
relation between stellar mass and dark matter halo mass in 
the models. The scatter around the mean gas fraction at 
fixed dark matter halo mass is set by a) the spin parame- 
ter of the halo, which determines the contraction factor of 
the accreted gas and hence the rate at which it will be con- 
sumed into stars, and b) the recent gas accretion history of 
the galaxy. The fact that the scatter in the models and data 
agree so well lends support to this basic picture. 

We note that in recent years, there have been a number 
of papers questioning the treatment of disk formation in 
semi-analytic models (e.g. Sales et al 2009). The primary 
objection is the assumption that the angular momentum of 
a galaxy expressed in units of that of its surrounding halo 
{jd = Jgai/ Jvir) , Correlates with mass of the galaxy divided 
by the mass of its dark matter halo (nid ~ Mgai/M^ir) in 
a manner that is insensitive to feedback, i.e. jd = nid is 
assumed for all halos. In our scheme this is not the case. 
We only assume that jd ~ nid for the gas that is cooling 
instantaneously at any given time in any given halo. Because 
we track the growth of galaxies and dark matter halos with 
time, and because feedback effects will affect the mass of gas 
that cools in halos of different masses at different redshifts, 
the resulting relation between the angular momentum of the 
galaxy and its mass fraction will depend on feedback. 

Figure 5 also shows that the models miss a minority 
population of galaxies with high atomic gas fractions and 
with low stellar masses, surface mass densities and concen- 
trations. What are these missing galaxies? In three recent 
papers, we have studied the properties of galaxies in our sur- 
vey that contain significantly more atomic gas than would 
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Figure 4. Scatter plots of H2 mass fraction (log M(H2)/M» ,left) 
and HI mass fraction (log M{HI)/M*, right) for COLD GASS 
galaxies are shown as blue points. Results from the FIO mod- 
els with KMT atomic-to-molecular gas prescription are shown as 
black points. 



be predicted from their UV/optical colours and sizes (Moran 
et al 2010; Wang et al 2011; Moran et al 2011). The main 
conclusion is that such galaxies have very blue outer disks 
with young stellar populations and low gas-phase metallici- 
ties. These results led us to propose that these unusually HI- 
rich systems have experienced a recent gas accretion event, 
resulting in growth of the outer disk. 

In the FIO models, the accreted gas is always added to 
the disk with an exponential surface density profile. In real 
galaxies, the accreted gas may initially be distributed in the 
outer regions of the galaxy. Dynamical perturbations in the 
form of spiral density waves, interaction with companions 
or with the non-axisymmetric gravitational potential of the 
surrounding dark matter halo, will eventually the cause the 
gas to flow inwards and reach high enough densities to form 
molecular clouds. If we change the gas accretion prescrip- 
tions so that accreted gas is initially dumped in the out- 
skirts of the disk, this may produce a tail of Hl-rich galax- 
ies. Observations of SFR, metallicity and gas profiles will be 
required to constrain models of this nature. 

3.2 The nature of the quenched population 

In this section, we investigate whether the "transition" be- 
tween the population of galaxies with detectable atomic and 
molecular gas and the populations with gas fractions be- 
low ~ 0.015 occurs in the same way in models as in the 
data. In Figures 7 and 8, we plot the fraction of HI/CO- 
detected galaxies as a function of stellar mass, stellar mass 
surface density, concentration index and specific star forma- 
tion rate. Results from the survey are shown as solid blue 
lines. The dashed blue lines indicate the zbla uncertainty 
in the detected fraction as a function of these parameters. 
Models results are shown in red and black (the colour-coding 
is the same as in Figures 2 and 3). The main conclusion 
from these two figures is that the models and the data do 
not agree. First, the dependence of the detected fraction on 
stellar mass is weaker in the data than in the models. This is 
true for both the HI and the CO-detected fractions. Second, 
the fraction of CO-detected galaxies drops very strongly as 
a function of both stellar surface density and concentration 
in the real data. This is not seen in the models. 

It is also instructive to look at the distribution of de- 
tected and non-detected galaxies in two-dimensional planes 
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Figure 5. Distribution functions of HI mass fraction (log 
M(HI)/M*) for COLD GASS galaxies with HI line detections in 3 
ranges of stellar mass, stellar mass surface density and concentra- 
tion (blue histograms). Error bars are calculated using bootstrap 
resampling. Results for the FIO models are shown as black and 
red curves. The colour-coding of the model curves has the same 
meaning as in Figures 2 and 3. 
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Figure 6. Distribution functions of H2 mass fraction (log 
M(H2)/M,) for COLD GASS galaxies with CO line detections 
in 3 ranges of stellar mass, stellar mass surface density and con- 
centration (blue histograms) . Error bars are calculated using 
bootstrap resampling. Results for the FIO models are shown as 
black and red curves. The colour-coding of the model curves has 
the same meaning as in Figures 2 and 3. 
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Figure 7. The fraction of COLD GASS galaxies with HI line 
detections is plotted as a function of stellar mass, stellar mass sur- 
face density, concentration and specific star formation rate (blue 
curve). Poisson errors are shown as blue dashed lines. We note 
that these are computed for independent bins, so the Icr errors 
can be read off directly. Results for the FIO models are shown as 
black (KMT atomic-to-molecular gas prescription) and red (BR 
pressure prescription) curves. 
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Figure 8. The fraction of COLD GASS galaxies with CO line 
detections is plotted as a function of stellar mass, stellar mass sur- 
face density, concentration and specific star formation rate (blue 
curve). Poisson errors are shown as blue dashed lines. Results for 
the FIO models are shown as black (KMT atomic-to-molecular 
gas prescription) and red (BR, pressure prescription) curves. 



of stellar mass and galaxy structural parameters. This is 
illustrated in Figures 9 and 10. We plot HI/CO detected 
galaxies in blue and HI/CO non-detected galaxies in red in 
the stellar surface density versus stellar mass, concentration 
versus stellar mass, and concentration versus stellar surface 
density planes. Results from the survey are shown in the 
top panels and results for "mock catalogues" of the same 
size generated from the simulations are shown in the bot- 
tom panels. In the data, there are clear thresholds in both 
(~ 3 X 10* M© kpc"^) and C (~ 2.6), that demarcate 
the location of almost all galaxies without detectable gas. 
No such threshold is seen in stellar mass M, . This is true 
for both the HI and the CO non-detections. 

In the models, the parameters that most clearly demar- 
cates the location of the HI non-detections are stellar mass 
surface density and stellar mass. It is the behaviour of the 
CO non-detections, however, that is most discrepant with 
the models. Galaxies without H2 are not confined to a spe- 
cific location in structural parameter space in the same way 
as in the data. 

We also note that the tight correlation between stel- 
lar surface density and concentration index seen in the top 
right panels of both figures, is not present in the models. 
In addition, the distribution of concentration indices in the 
models extends to much lower values than in the real data. 
The tight correlation between C and log ^* tells us that in 
the real Universe, galaxies with larger bulge-to-disk ratios 
have disks with higher stellar surface densities. 

In the FIO models, the effects of gas inflows on the disk 
are not taken into account, which may explain why there is 
no clear relation between bulge-to-disk ratio and stellar sur- 
face density. In a scenario where bulges form when gas flows 
inwards as a result of bar-driven inflows or other dynamical 
instabilities, the gas flows will not only form the bulge, but 
also increase the stellar surface mass density in the inner 
disk. 



3.2.1 Origin of quenching thresholds in the models 

We now elucidate the origin of the strong trend in the frac- 
tion of quenched galaxies as a function of stellar mass that is 
seen for the model galaxies. In current semi-analytic models, 
there are two physical processes that remove the supply of 
new gas to a galaxy and shut down star formation. 

(i) Radio Mode Feedback. Black holes are able 
to grow by accreting hot gas from the surround- 
ing halo. The growth rate is given by Mbh ~ 
K(/hot/0.1)(K»r/200km/s)3(MBH/10«MQ)MQ/yr (Cro- 
ton et al 2006; Guo et al 2011), where fhot is the ratio of 
hot gas mass to dark matter mass in the surrounding halo 
or sub-halo, Vvir is the virial velocity of the halo, Mbh is 
the black hole mass, and k is an efficiency parameter. Some 
fraction of the rest mass energy of the accreted material 
is assumed to be transferred to the surrounding hot gas 
by radio jets. The models assume an energy input rate 
£^radio = O.IMbhc^, whcrc c is the speed of light. This leads 
to a reduction in the cooling rate of hot gas of the form 

Afcool,c£F ~ Afcool — 2i?radio/^ir- 

(ii) Gas stripping. When a galaxy is accreted by a more 
massive dark matter halo, it becomes a "satellite". In the 
models of De Lucia & Blaizot (2007), the hot gas surround- 
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Figure 9. Top: COLD GASS galaxies are plotted in the 2-dimensional planes of stellar surface density versus stellar mass, concentration 
versus stellar mass and concentration versus stellar surface density. Galaxies with HI line detections are plotted in blue, while those without 
HI line detections are plotted in red. Bottom: Detections and non-detections from from our FIO model "mock catalogues" are plotted in 
the same 3 planes. 
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Figure 10. Top: COLD GASS galaxies are plotted in the 2-dimensional planes of stellar surface density versus stellar mass, concentration 
versus stellar mass and concentration versus stellar surface density. Galaxies with CO line detections are plotted in blue, while those 
without CO line detections are plotted in red. Bottom: Detections and non-detections from from our FIO model "mock catalogues" are 
plotted in the same 3 planes. 
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ing the satellite is stripped instantaneously, leading to a 
sharp reduction in the cooling rate onto the satellite. In the 
more recent models of Guo et al (2011), the dark matter and 
hot gas surrounding the satellite are removed more gradually 
both by tidal forces and by ram-pressure stripping. 

The parameter k is tuned to reproduce the high mass 
end of the stellar mass function. The net effect is illustrated 
in Figure 11 where we plot galaxies from our model cata- 
logue in the 2D planes of subhalo mass versus log M, , sub- 
halo mass versus log /x* , and subhalo mass versus C. Galax- 
ies that are predicted to be detected in CO are colour-coded 
blue, while those predicted to be non-detections are colour- 
coded red. Figure 11 shows a clear transition between active 
and quenched galaxies at a subhalo mass of lO^^Afo, which 
is a factor of ~ 3 larger than the mass where dark matter 
halos are predicted to transition to hosting a static halo of 
hot gas (Birnboim & Dekel 2003; Dekel & Birnboim 2006; 
Croton et al 2006). The reason why the transition is very 
sharp in halo mass is because the primary dependence of 
the mass that is accreted by the black hole is on halo virial 
velocity and not on black hole mass. We note that there is 
an additional population of quenched galaxies in low mass 
subhalos. These are the satellite galaxies, where the subhalo 
has been stripped by tidal forces. As can be seen, they form 
a small minority of the quenched population in the stellar 
mass range of the galaxies in the GASS survey. 

Figure 11 clearly illustrates why quenching in the model 
galaxies is most strongly dependent on stellar mass and sur- 
face density, and largely independent of concentration. This 
is because in the model, stellar mass and surface density cor- 
relate with the virial mass of the host halo, but concentration 
does not. There is direct observational evidence in support 
of these results. Mandelbaum et al (2006) used weak gravita- 
tional lensing to derive halo masses of galaxies as a function 
of stellar mass and morphology. They found stellar mass to 
be a good proxy for halo mass. For a given halo mass, the 
stellar mass was independent of concentration/morphology 
below M* = 10" Mq. 

3.2.2 Origin of quenching thresholds in real galaxies 

In the data, quenching thresholds are seen as a function 
of stellar surface density and concentration, but not as a 
function of stellar mass. Even for galaxies with stellar masses 
as high as 10^^ Mq, which are observed to reside in halos 
with masses in the range 3 x lO^'^ — 1O^^M0, HI and CO are 
generally still detected if concentrations and densities are 
low. This strongly suggests that processes associated with 
bulge formation must be responsible for shutting off the gas 
supply in galaxies. 

We note that this hypothesis was already put forward 
by Kauffmann et al (2006), based on an analysis of the scat- 
ter in the colours and spectral properties of galaxies as a 
function of stellar mass, stellar surface density and concen- 
tration. Kauffmann et al (2006) made an ansatz that scatter 
in specific star formation rate reflected scatter in gas con- 
tent. The fact that specific star formation rates as well as 
their scatter decreased sharply above a characteristic density 
threshold of 3 x 10* Mq kpc~'^ and concentration index of 2.6 
and that this threshold was largely independent of the stellar 
mass of the galaxy, was taken as evidence that gas accretion 




D(4000) 



Figure 12. The Lick H(5a index is plotted as a function of 
4000 A break strength for COLD GASS with CO line detections 
and log fit > 8.6 (blue). CO line detections and log^* < 8.6 
(cyan) and for COLD GASS galaxies with no CO line detections 
(red). The errors on each measurement are also indicated for each 
galaxy. 

was no longer occurring in bulge-dominated systems. The 
COLD GASS and GASS surveys have now demonstrated 
that the neutral gas content of galaxies is weakly dependent 
on stellar mass and decreases sharply at ^, > 3 x 10* 
kpc"^ and C > 2.6. 

We have not yet ascertained why bulge-dominated 
galaxies no longer accrete gas and form stars efflciently. Fig- 
ure 3 shows that galaxies with /x* > 3 x 10* Mq kpc~'^ have 
molecular gas mass fractions that are slightly depressed rel- 
ative to model predictions. Could such galaxies be transi- 
tioning to the red sequence on short timescales as proposed 
by Schawinski et al (2009)? 

We note that molecular gas is generally concentrated 
towards the inner regions of galaxies. The SDSS fiber spec- 
tra probe the central 1-2 kpc central regions of the GASS 
galaxies and contain diagnostics of the recent star formation 
history in the bulge. In Figure 12, we plot H(5a as a function 
of D„(4000) for the galaxies in the survey. Galaxies where 
CO was not detected are colour coded red, those with de- 
tections and log fj., > 8.6 are colour coded-blue, and those 
with detections and log /x* < 8.6 are colour-coded cyan. As 
discussed in detail in Kauffmann et al (2003a) , by combining 
these two stellar absorption line indices, we can diagnose if 
the recent star formation histories of a population of galax- 
ies have been smooth or "bursty" on average. Starbursts of 
duration less than 1-2 hundred million years, displace a sig- 
nificant fraction of galaxies to higher values of H5a at a given 
value of D„(4000). Likewise, if star formation is truncated 
over a timescale of less than a few hundred million years, 
galaxies will also be displaced to higher than average values 
of H5yi for around a Gyr following the truncation event (see 
for example, Kauffmann et al 2004). 

As can be seen, most GASS galaxies lie on the same 
locus in the plane of US a versus D„(4000). There is no ev- 
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Figure 11. FIO model galaxies are plotted in the 2-dimensional planes of subhalo mass versus stellar mass, subhalo mass versus stellar 
surface density and subhalo mass versus concentration. Galaxies predicted to have CO line detections are plotted in blue, while those 
predicted not to have CO line detections are plotted in red. 



idence of any net displacement as a function of stellar sur- 
face density or as a function of H2 mass fraction, indicating 
that the star formation histories of most high surface den- 
sity galaxies with gas have been smooth. We note, however, 
that galaxies with CO line detections and high stellar sur- 
face densities do span the largest range in stellar population 
parameters. The have D„(4000) values as low as those of 
the low-density H2-rich population, and as high as those 
of the "quenched" population. We also note that the three 
post-starburst galaxies that are clearly displaced to higher 
valued of H5a, are all detected in H2 have have high stellar 
surface densities. All of them have unusually high molecular- 
to-atomic gas ratios, and one is clearly the end product of a 
recent merger (tidal features are visible in the SDSS image). 

This region of parameter space is clearly quite com- 
plex, and may consist of several distinct sub-populations of 
galaxies on different evolutionary trajectories. If so, addi- 
tional data will be required to pull these apart. In particu- 
lar, studies of the distribution and the kinematics of the gas 
may shed significant light on its origin and its eventual fate. 



4 SUMMARY 

We compare the semi-analytic models of galaxy formation 
of Fu et al. (2010), which track the evolution of atomic and 
molecular gas in galaxies, with gas fraction scaling relations 
derived from a stellar mass-limited sample of 299 galaxies 
from the COLD CASS survey. These galaxies have mea- 
surements of the CO (1-0) line from the IRAM 30-m tele- 
scope and the HI line from Arecibo, as well as measurements 
of stellar masses, structural parameters and star formation 
rates derived from GALEX-I-SDSS photometry. 

Our analysis addresses two key questions: 1) Can the 
semi-analytic disk formation models explain the observed 
gas fraction scaling relations in galaxies with gas and on- 
going star formation?, 2) Does the transition between the 
population of galaxies with gas and the "quenched" popu- 
lation without gas occur in the same way in the models and 
in the data? 

In answer to the first question, we conclude that the 
disk models provide a reasonable description of how con- 
densed baryons are partitioned into stars, atomic gas and 
molecular gas as a function of galaxy mass and size. Trends 



as a function of bulge-to-disk ratio are not well reproduced. 
In particular, the models do not account for the fact that at 
fixed stellar mass, there is a tight relation between the size 
of a galaxy and its bulge-to-disk ratio. 

In answer to the second question, we conclude that 
our data disagree with the current implementation of radio- 
mode feedback in the models. In the models, the observ- 
able parameter that best predicts whether a galaxy has been 
quenched is its stellar mass. Our data shows that the fraction 
of quenched galaxies is largely independent of stellar mass, 
but depends strongly on galaxy bulge-to-disk ratio and stel- 
lar surface density. In other words, even low mass galaxies 
with bulges have high probability of being quenched. We 
conclude that processes associated with bulge formation are 
thus likely to be responsible both for depleting the neutral 
gas in galaxies and for shutting off the gas supply in these 
systems. 

Solving these problems will require substantive changes 
to the way gas transport and bulge and black hole formation 
is treated in the models, as well as to the way feedback from 
AGN is implemented. We now outline our best guess as to 
how this might work in practice. 

In the current models, newly accreted gas is assumed 
to have an exponential profile. If the profile is initially shal- 
lower than exponential, more atomic gas would collect in an 
largely inert reservoir in the outer regions of the galaxy. This 
gas would later be driven towards the center of the galaxy 
by dynamical perturbations, where it would form molecular 
gas and stars in the inner region of the disk. The inflow- 
ing gas would also contribute to the growth of the bulge. A 
prescription of this kind may result in a better correlation 
between stellar surface mass density and bulge-to-disk ratio. 

There is now evidence from studies of complete sam- 
ples of nearby galaxies that bars are very common in disk 
galaxies. Barazza, Jogee & Marinova (2008) find that 70% of 
disk-dominated galaxies host bars. Moreover, barred galax- 
ies (in particular galaxies with strong bars) are associated 
with higher central molecular gas fractions and enhanced 
rates of central star formation (Sakamoto et al 1999; Ellison 
et al 2011; Wang et al 2012, in preparation), suggesting that 
bulges are forming in these systems . Finally, barred galaxies 
do not exhibit any axcess of nearby companions, suggesting 
that they are not triggered by interactions (Li et al. 2009). 

All of this suggests that internally-driven gas transport 
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processes are important in the formation of bulges, but re- 
solved maps of the atomic and molecular gas distributions 
in complete samples of galaxies and more detailed theoreti- 
cal calculation of gas inflow rates in galaxies are required in 
order to quantify this in more detail. 

In the current models, radio-mode feedback suppresses 
the cooling of gas in halos more massive than 10^^ Mq. In 
addition, in massive, high density galaxies, atomic gas is 
transformed very efficiently into molecular gas and thence 
into stars. These two mechanisms appear to be sufficient to 
place all "Hl-quenched" galaxies in the high M« amd p» 
corner of parameter space, as seen in the bottom panel of 
Figure 9. The same mechanisms are not, however, sufficient 
to place the "H2-quenched" galaxies in this same region of 
parameter space. Nor can they explain why galaxies deficient 
in both HI and H2 are almost always found in galaxies with 
substantial bulge components. The simplest inference is that 
the neutral gas in galaxies is either consumed or removed 
when bulges form, and that subsequent gas accretion into 
the atomic phase is suppressed in a significant fraction of 
bulge-dominated galaxies. 

More detailed studies of the spatially-resolved kinemat- 
ics of the gas and how gas motions across galaxies correlate 
with the presence of accreting black holes, jets and the dis- 
tribution of star formation, will be needed before we are able 
to pinpoint the actual physical mechanisms responsible for 
removing/depleting the gas in these systems (see for exam- 
ple Hopkins et al 2011). Quenching may occur over short 
timescales during particular phases of bulge formation, so it 
may be necessary to survey large samples before a definitive 
conclusion can be reached. In addition, the mechanisms that 
inhibit the late accretion of gas in bulge-dominated systems 
remain to be clarified. 
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